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Abstract 



We are given a stack of pancakes of different sizes and the only allowed operation 
is to take several pancakes from top and flip them. The unburnt version requires 
the pancakes to be sorted by their sizes at the end, while in the burnt version they 
additionally need to be oriented burnt-side down. We present an algorithm with 
the average number of flips, needed to sort a stack of n burnt pancakes, equal to 
7n/4 + 0(1) and a randomized algorithm for the unburnt version with at most 
17n/12 + 0(1) flips on average. 

In addition, we show that in the burnt version, the average number of flips of 
any algorithm is at least n + 0(n/logn) and conjecture that some algorithm can 
reach n + G(n/ log n). 

We also slightly increase the lower bound on g(n), the minimum number of flips 
needed to sort the worst stack of n burnt pancakes. This bound together with the 
upper bound found by Heydari and Sudborough in 1997 gives the exact number 
of flips to sort the previously conjectured worst stack —I n for n = 3 (mod 4) and 
n > 15. 

Finally we present exact values of f{n) up to n = 19 and of g{n) up to n = 17 
and disprove a conjecture of Cohen and Blum by showing that the burnt stack — /15 
is not the worst one for n = 15. 
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1 Introduction 



The pancake problem was first posed in [I]. We are given a stack of pancakes each two 
of which have different sizes and our aim is to sort them in as few operations as possible 
to obtain a stack of pancakes with sizes increasing from top to bottom. The only allowed 
sorting operation is a " spatula flip" , in which a spatula is inserted beneath an arbitrary 
pancake, all pancakes above the spatula are lifted and replaced in reverse order. 

We can see the stack as a permutation 7T. A flip is then a prefix reversal of the 
permutation. The set of all permutations on n elements is denoted by S n , f(n) is the 
minimum number of flips needed to obtain (1, 2,3, ... ,n) from 7r and 

f(n) := max f(7i). 

The exact values of f(n) are known for all n < 19, see Table [T] for their list and 
references. In general 15 Lrz/14j < f(n) < 18n/ll + 0(1). The upper bound is due to 
Chitturi et al. [2] and the lower bound was proved by Heydari and Sudborough [U] . These 
bounds improved the previous bounds 17n/16 < f{n) < (5n + 5)/3 due to Gates and 
Papadimitriou [6J, where the upper bound was also independently found by Gyori and 
Turan [7]. 

A related problem in which the reversals are not restricted to intervals containing the 
first element received considerable attention in computational biology; see e. g. [8]. 

A variation on the pancake problem is the burnt pancake problem in which pancakes 
are burnt on one of their sides. This time, the aim is not only to sort them by their sizes, 
but we also require that at the end, they all have their burnt sides down. Let C = (tt, v) 
denote a stack of n burnt pancakes, where 7r 6 S n is the permutation of the pancakes 
and v G {0, 1}" is the vector of their orientations (f« = if the i-th pancake from top is 
oriented burnt side down). Pancake i will be represented by i if its burnt side is down 
and i if up. Let 

/l\ /I\ 

and — /„ = 

\nj \n J 

Let g{C) be the minimum number of flips needed to obtain I n from C and let 



g{n) : = max g{{ir,v)). 

TT€S n ,v€{0,l} n 



Exact values of g(n) are known for all n < 17, see Table [TJ In 1979 Gates and 
Papadimitriou [H] provided the bounds 3n/2 — 1 < g(n) < 2n + 3. Since then these were 
improved only slightly by Cohen and Blum [3] to 3n/2 < g(n) < 2n — 2, where the upper 
bound holds for n > 10. The result g(lQ) = 26 further improves the upper bound to 
2n — 6 for n > 16. Cohen and Blum also conjectured that the maximum number of flips 
is always achieved for the stack — I n . But we present two counterexamples with n = 15 
in Section [6j 
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The stack — I n can be sorted in (3(n + 1))/2 flips for n = 3 (mod 4) and n > 23 [5]. In 
Section [3] we present a new formula for determining a lower bound on the number of flips 
needed to sort a given stack of burnt pancakes. The highest value that this formula gives 
for a stack of n pancakes, is |_(3(n+l))/2j for the stack — I n . These bounds together with 
the known values of g(—I 15 ) and g(— Jig) give g(—I n ) = (3(n + l))/2 if n = 3 (mod 4) 
and n > 15. 
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n = 3(mod4) L^J Corollary 1 

Table 1: known values of f(n), g{n) and g(—I n ) 

We present an algorithm that needs on average 7n/4 + 0(1) flips to sort a stack 
of n burnt pancakes and a randomized algorithm for sorting n unburnt pancakes with 
Yin/ 12 + 0(1) flips on average. We also show that any algorithm for the unburnt version 
requires on average at least n — 0(1) flips and in the burnt version n + Q(n/ logn) flips 
are needed on average. Section [7] introduces a conjecture that the average number of flips 
of the optimal algorithm for sorting burnt pancakes is n + G(n/logn). 



2 Terminology and notation 

The stack obtained by flipping the whole stack C is C. The stack — C is obtained from 
C by changing the orientation of each pancake while keeping the order of pancakes. 
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If two unburnt pancakes of consecutive sizes are located next to each other, they are 
adjacent. Two burnt pancakes located next to each other are adjacent if they form a 
substack of I n or of /„. Two burnt pancakes located next to each other are anti-adjacent 
if they form a substack of — I n or of —I n . 

In both versions a block in a stack C is an inclusion-wise maximal substack S of C 
such that each two pancakes of S on consecutive positions are adjacent. A substack S 
of a stack C with burnt pancakes is called a clan, if —S is a block in —C. Pancake not 
taking part in a block or a clan is free. 

If the top i pancakes are flipped, the flip is an i-flip. 



3 Lower bound in the burnt version 

Theorem 1. For each n 



g(-in) > 



3 71+1 



Proof. 

The claim is easy to verify for n < 2, so we can assume n > 3. 

A block (clan) is called a surface block (clan) if the topmost pancake is part of it, 
otherwise it is deep. 

We will assign to each stack C the value v(C): 

v(C) :=a(C)-a-(C)- l -(b(C)-b-(C)) + ^ 
where 

a(C) := number of adjacencies 

b(C) := number of deep blocks 

1 if the pancake on top of the stack is the free 1 or 
o(C) := < if 1 is in a block (necessarily with 2) 
otherwise 



1(C) := 

11(C) :-- 

a-(C) 
b-(C) 

0- (C) 

1- (C) 
11- {C) 



J 1 if the lowest pancake is n 
\ otherwise 

J 1 if the lowest pancake is n and the second lowest is n — 1 
\ otherwise 

a(—C) = number of anti-adjacencies in C 

b(—C) = number of deep clans in C 

o(-C) 

l(-C) 

ll(-C). 
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Lemma 2. If C and C are stacks of at least two pancakes and C can be obtained from 
C by a single flip, then 

Av := v(C) - v(C) < \. 
Therefore the minimum number of flips needed to sort a stack C is at least 

-Ml n ) ~ V{C)) 



Proof. 

First we introduce notation for contributions of each of the functions to Av: 



Aa 


:= a(C) - a(C) 


Aa 


Ab 


:=- l -(b(C')-b(C)) 


Ab 


Ao 


:=\{o{C') -o(C)) 


Ao 


Al 


:= 1(C) - 1(C) 


Al 


All 


■■=1(11(0')- 11(C)) 


All 



: (b-(C')-b-(C)) 



1 



o-(C')-o-(C)) 



Ai- ■= -(r(C)-r(c)) 
= -\(ir(C)-ir(C)) 



Observation 3. Values of Aa, Aa , Al and Al are among {0,1,-1}. Values of Ab, 
Ab~ , Ao, Ao~ ; All and All~ are among {0,1/3,-1/3}. 

Proof. The only nontrivial part is Ab < 1/3 and symmetrically Ab~ < 1/3. For contra- 
diction suppose Ab > 1/3, which can only happen when one block was split to two free 
pancakes and another block became surface in a single flip. But the higher of the two 
pancakes that formed the split block will end on top of the stack after the flip. Therefore 
no block became surface. To show Ab~ < 1/3 we consider the flip : — C — > — C, for 
which 

I > A,6 = -\(b(-C) - b(-C')) = -\(b~(C) - b-(C')) = \(b-(C>) - b-(C)) = Ab-. 



□ 



The proof of the lemma is based on restricting possible combinations of values of the 
above defined functions. 

• Both Al and Al~ are positive. This would require the pancake n to be before and 
after the flip at the bottom of the stack each time with a different orientation. But 
this is not possible when n > 1. 
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• Exactly one of Al and Al is positive. The case Al > can be transformed to the 
case Al > by considering the flip <fi : — C — > — C , for which 



A^v 


:=v{ 


-c)- 


-v(-C') = -v(C) - (-v(C)) = v(C) 


-v(C) = Av, 


A^l 


:=/(- 


-c)- 


i(-c') = r(c) - r(c") = -(r(C) 


-r(c)) = a/- 


A^r 


:= r 


(-c) 


- r(-c') = al 





The equality v(—C) = —v(C) follows from the definition of v(C). 

If the value of / changes, the flip must be an n-flip. Therefore Aa = Aa~ = 0. 
Because Al — 1, the pancake n has to be at the bottom of the stack after the flip, 
so All~ = 0. Moreover neither a clan nor the pancake 1 could be on top of the stack 
before the flip so Ab~ < and Ao~ < 0. Because All =1/3 implies a block on top 
of the stack before the flip and Ao =1/3 implies no block on top of the stack after 
the flip, we obtain 

All = - & Ao < =>- Ab < 0, 

All < & Ao = i Ab < 0, 

All = - LAo=-^Ab< --. 
3 3 ~ 3 

In any of the cases All + Ao + Ab < 1/3 and Av < 4/3. 
From now on, we can assume Al, AZ~ < 0. 

• At least one of All and All~ is positive. If both of them were positive then again the 
pancake n would be at the bottom of the stack before and after the flip, each time 
with a different orientation. Similarly to the previous case, we can choose All~ = 
and All = 1/3. Because Al < 0, the last flip was an (n — l)-flip, the pancake at 
the bottom of the stack is n and the pancake on top of the stack before the flip was 
(n- 1). Therefore Aa = 1, Aa~ = 0, Ao~ < and Ab~ < 0. 

If pancake n — 1 was part of a block before the flip, then this block became deep, 
otherwise pancakes n — 1 and n created a new deep block. Thus Ab < 0. No block 
was destroyed and if Ao = 1/3, then no block became surface and thus Ab = —1/3. 
All in all Av < 4/3. 

In the remaining cases we have Al, A/ - , All, All~ < 0. 

• Both Ao and Ao - are positive. Because Ao - > then either 1 was in a clan or on 
top of the stack with burnt side down before the flip. If 1 was in a clan, then a single 
flip would not make it either a part of a block or a free 1 on top of the stack and 
thus Ao would not be positive. Using a similar reasoning for Ao, we obtain that the 
flip was a 1-flip, the topmost pancake before the flip was 1 and the second pancake 
from top is different from 2. Thus Aa = Aa~ = Ab = Ab~ = and Av < 2/3. 
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Exactly one of Ao and Ao is positive; without loss of generality it is Ao. This can 
happen only in two ways. 

— We did an i-flip, the topmost pancake before the flip was 2 and the (i + l)-st 
pancake is 1. Then Aa = 1, Aa~ = 0, Ab < and Ab~ < and so Av < 4/3. 

— We did an i-flip, the i-th pancake before the flip was 1 and neither the (i — l)-st 
nor the (i + l)-st pancake was 2. Then Ab < and Aa - < 0. If Aa < 0, then 
Av < 2/3, otherwise Ao" < and Av < 4/3. 

Now only Aa, Aa~, Ab and Ab~ can be positive. 
If Aa = Aa - = 1, then the flip was either 



z + 1 



V 



V 



or 



/ 



i - 1 

i 



V 



/ 



l + l 

i 



V 



/ 



In both cases the topmost pancake before the flip was not part of a clan and the 
topmost pancake after the flip is not part of a block, so the number of deep blocks 
increased and the number of deep clans decreased and Aw < 4/3. 

Exactly one of Aa and Aa~ is positive; without loss of generality Aa = 1, Aa~ < 0. 
Neither a new clan was created, nor became deep, so Ab~ < and Av < 4/3. 



None of Aa and Aa is positive, so Av < 2/3. 



n + 2/3 and v(—I n 
number of flips needed to transform — I n to /„ is at least 

"3 



It is easy to compute that v(I n 
os needed to tran 

(v(I n )-v(-Q) 



-n 



□ 

2/3 and thus the 









"3 




3(n+ 1) 
2 




f!M)l 




-n + 1 





□ 

Corollary 4. For all integers n > 15 with n = 3 (mod 4), 

3(n+i; 

9{-In) = g 

Proof. The lower bound comes from Theorem [TJ For all n > 23 with n = 3 (mod 4), the 
upper bound was proved by Heydari and Sudborough [9J. The exact value for n — 15 
was computed by Cohen and Blum [3] and the exact value for n = 19 is computed in 
Section E □ 
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4 Algorithm for the burnt version 



In this section we will design an algorithm that sorts burnt pancakes with small average 
number of flips. 

First we will show a lower bound on the average number of flips of any algorithm that 
sorts a stack of n burnt pancakes. 

Theorem 5. Let av opt (n) be the average number of flips of the optimal algorithm for 
sorting a stack of n burnt pancakes. For any n > 16 

n 3 

av opt {n) > n + — -. 

16 log 2 n 2 

Proof. We will first count the expected number of adjacencies in a stack of n burnt 
pancakes. A stack has n — 1 pairs of pancakes on consecutive positions. For each such 
pair of pancakes, there are 4n(n — 1) equally probable combinations of their values and 
orientations and the pancakes form an adjacency in exactly 2(n — 1) of them. From the 
linearity of expectation 

, , / . 1 In - 1 
E M = („-!)- = -_. 

Therefore at least half of the stacks have no adjacency. 



• First we take a half of the stacks such, that it contains all the stacks which have 
some adjacency. The stacks of this half have less than 1 adjacency on average. Each 
flip creates at most one adjacency, therefore when we want to obtain the stack I n 
with n — 1 adjacencies, we need at least n — 2 flips on average. 

• The other half contains n\ ■ 2 n_1 stacks each with no adjacency, thus requiring at 
least n — 1 flips. For each stack we take one of the shortest sequences of flips 
that create the stack from I n and call it the creating sequence of the stack. Note 
that creating sequences of two different stacks are different. We will now count the 
number of different creating sequences of length at most n — 1 + n/(4 log 2 n), which 
will give an upper bound on the number of stacks with no adjacency that can be 
sorted in n — 1 + n/(41og 2 n) flips. Shorter creating sequences will be followed by 
several 0-flips, therefore we will consider n + 1 possible flips. A split-flip is a flip in a 
creating sequence that decreases the number of adjacencies to a value smaller than 
the lowest value obtained before the flip. Therefore there are exactly n — 1 split-flips 
in each of our creating sequences. In a creating sequence, the i-th split-flip removes 
one of n — i existing adjacencies and therefore there are n — i possibilities how to 
make the i-th split-flip. The number of different creating sequences of the above 
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given length is at most 
'n - 1 + 



11 



log 2" • (n-1)!- ( n + l) n/(41og2n) 



41og 2 n 

(\ n/(4 log 2 n) 
+ — - — ) " ■ (n- 1)! ■ (2n) n/(41og2n) 

< (n - 1)! • (2n) n/(41og2n) • (2n) n/(41og2n) 
<(n-l)!-(n 5 / 4 ) 2n/(41 ° S2n) 

< (n- 1)! -2 5n/8 

< ^! -2 n . 
4 

Thus at least half of the stacks with no adjacency need more than n—l+n/ (4 log 2 n) 
flips while the rest needs at least n— 1 flips. Therefore in this case the average number 
of flips is at least 

n 

n - 1 + 



8 log 2 n 

The overall average number of flips is then 

i \ 3 n 

avoptin) >n-- + 



2 16 log 2 n 

□ 

Theorem 6. There exists an algorithm that sorts a stack of n burnt pancakes with the 
average number of flips at most 

7 

-n + 5. 
4 

Proof. Let C n denote the set of all stacks of n burnt pancakes, h(C) will be the number 
of flips used by the algorithm to sort the stack C and let 



H(n) := h (C), 

cec„ 

H(n) H{n) 



av[n 



|C n | 2n\C, 



n-l 



The algorithm will never break previously created adjacencies. This allows us to con- 
sider the adjacent pancakes as a single burnt pancake. In each iteration of the algorithm 
one adjacency is created, the two adjacent pancakes are contracted and the size of the 
stack decreases by one. We stop when the number of pancakes is two and the algorithm 
can transform the stack to the stack (1) in at most four flips. 

However for the simplicity of the discussion, we will not do such a contraction for 
adjacencies already existing in the input stack (as can be seen in the proof of Theorem [5], 
there are very few such adjacencies, so the benefit would be negligible). 
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One more simplification is used. Before each iteration, the algorithm looks at the 
topmost pancake and cyclically renumbers the pancakes so as to have the topmost pancake 
numbered 2 — pancake number j will become j + s + kn, where s = (2 — 7r(l)) and k 
is an integer chosen so as to have the result inside the interval {1, . . . , n}. Let be the 
set of stacks with n burnt pancakes and the pancake number 2 on top. When we end up 
with the stack (1), we in fact have 



i + l 



n 
1 
2 



Vizi!/ 

for some i G {1, 2, ... n}. This stack needs at most four more flips to become I n . Therefore 
civ (2) < 8. We will do four flips at the end even if they are not necessary. Then the 
number of flips will not be changed by a cyclic renumbering of pancakes and H{n) = 
n -Y.ceCi h (C)- 

• If the stack from can be flipped so that the topmost pancake will form an 
adjacency, we will do it: 



/ 2 
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X 
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1 




\ Y 
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V 


Y 


J 




i- 



or 



/ 2 \ / X \ / vi \ /vii 
X W 2 \ \ \ \_ 

3 3 \ Y' j \ Y" 

\yJ \Y J v / v 

Each stack from C n _i appears as a result of the above described process for exactly 
one stack from C^. 

• If no adjacency can be created in a single flip, we will look at both pancakes 1 and 
3 and analyze all possible cases. Note that this time when 2 has its burnt side up, 
then 3 has its burnt side up and similarly 2 implies 1. 
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7. 
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Again each stack from C n _i appears as a result of the process for exactly one stack 
from C^, but we needed two additional flips in two of the cases to ensure this. We 
did four flips in a quarter of the cases and two flips in all other cases. Each case has 
the same probability and hence the average number of flips is 5/2. 



All in all 



H(n)=n-( (MC) + 1)+ £ ( h ^ + D) 



= 2nH(n-l) + -ra|C n _i|, 



av{n) 



2nH(n-l) + fra|C„_i| 
2n|C n _i| 



au(n - 1) + - = av(2) + -(n - 2) < -n + 5. 



□ 



5 Randomized algorithm for the unburnt version 

Observation 7. Lei av' opt (n, 0) fre £/ie average number of flips of the optimal algorithm 
for sorting a stack of n unburnt pancakes. For any positive n 

av' opt {n,0) >n-2. 

Proof. We will now count the expected number of adjacencies in a stack of n pancakes. 
For the purpose of this proof we will consider the pancake number n at the bottom of 
the stack as an additional adjacency; this has probability 1/n. Pancakes on consecutive 
positions form an adjacency if their values differ by 1; the probability of this is 2/n. 
Therefore the expected number of adjacencies is 

1 2 
EM = - + (n- 1)- < 2. 

n n 
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Each flip creates at most one adjacency, therefore when we want to obtain the stack 
/„ with n adjacencies, the average number of flips is at least n — 2. □ 

Theorem 8. There exists a randomized algorithm that sorts a stack ofn unburnt pancakes 
with the average number of flips at most 

— n + 9, 
12 ' 

where the average is taken both over the stacks and the random bits. 

Proof. If two pancakes become adjacent, we contract them to a single burnt pancake; its 
burnt side will be the one where the pancake with higher number was. Therefore in the 
course of the algorithm, some of the pancakes will be burnt and some unburnt. For this 
reason we say that two pancakes are adjacent if the unburnt ones of them can be oriented 
so that the two resulting pancakes satisfy the definition of adjacency for burnt pancakes. 

Let V n> b denote the set of all stacks of n pancakes b of which are burnt and let b be 
the stacks from U nj ;, with the pancake number 2 on top. Let k(C) be the number of flips 
needed by the algorithm to sort the stack C and let 



K(n,b):= k (C), 



ceu„ >6 

av ( n ' fe ) := iiT i • 

When there are only two pancakes left, we can sort the stack in at most 4 flips. 
Similarly to the burnt version, we will sometimes cyclically renumber the pancakes. After 
renumbering them back at the end, we will do 4 flips to get the sorted stack. Therefore 
av'(l, 0) = av'CL, 1) = 4, av'(2, 6) < 8 for any b G {0, 1, 2} and K(n, b) =n-^ raj2 k(C). 

n,b 

The algorithm first cyclically renumbers the pancakes so as to have the topmost pan- 
cake numbered 2 thus obtaining a stack from b . Then we look at the topmost pancake. 
If it is unburnt, we uniformly at random select whether to look at 1 or 3; if it is burnt 
and the burnt side is down, we look at 1 and in the case when the burnt side is up, we 
look at 3. 

Notice that we could also look at both pancakes 1 and 3. But if we joined only two 
of the pancakes 1, 2 and 3 we would have to count the average number of flips for each 
combination not only of the number of pancakes and the number of burnt pancakes, but 
also of the number of pairs of pancakes of consecutive sizes exactly one of which is burnt. 
This would make the calculations too complicated. We could also join all three of them, 
but this would lead to a worse result. 

2 m 

I. Both the pancakes we looked at are unburnt. The set of such stacks is U n ^- Note 
that stacks with pancake 2 unburnt and exactly one of pancakes 1 and 3 unburnt 
belong to this set from 50% — with 50% probability, we choose to look at the 
unburnt pancake. Let av^ri, b) be the weighted average number of flips used by the 
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2 m 

algorithm to sort a stack from U n ~^, where the weight is the ratio with which the 
2 m 

stack belongs to U n ^- 



/ 2 \ 

X 
1 

\ Y J 



2 
1 

\Y ) 



X' 
I 
Y' 



G U 



n-1,6+1 



/ 2 \ 

X 
3 

V Y J 



2 
3 

V Y J 







f X" 








Y> / 




i Y" 



1,6+1 



For each stack from U n -i^+i there are exactly 6+1 its cyclic renumberings each 
appearing as a result with a 50% probability. Thus we can compute the average 
number of flips in this case: 

avj^n, b) = av'(n — 1, b + 1) + 1. 
II. The topmost pancake is unburnt, while the other pancake we looked at is burnt. 



/ 2 \ /X\ 
X 2 
T T 

V Y J \Y J 



X' 

^ I T 

Y' 



n— l.i 



( 2 \ / T_\ /X \ 

X X - 

1 ~" 2 2 

\Y J \Y J \Y J 



X' 

<=> | 1 | GU n _l, b 



The case when we looked at pancake 3 is similar, so we can conclude that 



avfrrln, b) = av'(n — 1,6) + -. 



III. The topmost pancake is burnt, while the other one we looked at is unburnt. 



/ 2 \ 

X 
3 

\Y J 



2 
3 

V Y J 







(X 








Y' 1 




V Y' 
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/ 2 \ fX\ 

X 2 
1 1 



Each stack from U n _i i b appears as a result exactly once for b its cyclic renumberings. 
Therefore 

cm^jj/n, b) = av'(n — 1,6) + 1. 



IV. Both the pancakes we looked at are burnt. In half of the cases the two pancakes can 
be joined in a single flip: 



/ 2 \ 

X 
3 

V Y ) 



2 
3 

\Yj 









2 


IH 


i 


Y> 1 







G U 



n-1,6-1 



/ 2 \ (X\ 

X 2 
T ~* T 



X" 

1 

Y" 



g U„_i. 



6-1 



Otherwise we need three flips to join the two pancakes: 



/ 2 \ 




( 1 \ 




( ^ ^ 






X 




X 




X 




3 


3 




3 




2 




2 


\Y ) 




\Y ) 















f X" 




IH 











G U n _i 6_i 



/ 2 \ 
1 

V ^ / 



/ 2 \ 

X 

1 

V Y J 



( }_\ 

X 

2 

v y y 



( x\ 
1 

2 

v y y 



X" 

1 

F" 



G U 



n-1,6-1 



a?JTT7|(n, 6) = av'(n — 1, 6 — 1) + 2. 



Altogether 



After summing up all the above average numbers of flips multiplied by their probabil- 
ities, we obtain: 
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For 1 < b < n 



i, ,n (n — b)(n — b—l) (n - 6)6 ( . . . . ,.\ 

GU (n ' } = n(n - 1) "W 1 ' } + n(n-l) W"' &) + ""EI]^' 6) J + 



6(6-1) 
n(ra — 1) 
(n-6)(n-6-l) 



+ — — — Q Wf n > fe ) 



n(n — 1) 



(l + au'(n- 1,6 + 1))+ 



+ 2^-^ f ^ + av'(n - 1, 6)) + ( 2 + ^'(« -1,6-1))- 

n(n-l)V 4 / ra(ra-l) 

• For 6 = 

Tl\Th — 1 ) 

ai/(n, 0) = —. ratlin, 0) = 1 + av'(n - 1, 1). 

n(n — 1) lAJ 

• For b = n 

av (n, n) = — rawprnn, n) = 2 + av (n — 1, n — 1). 

n(n — 1) LL^J 

Instead of solving these recurrent formulas, we will use them to bound av'(n, b) from 
above by the following function: 

+/ ^ 17 7 l l(n-6 + l)6 n 

at> + (n, 6) := — n H 6 — + 9. 

v ' ; 12 12 6 n 

Lemma 9. For any nonnegative n and b, such that b is not greater than n 

av + (n,b) > av'(n,b). 
Proof. We will use induction on the number of pancakes. 

• For n = lwe have av'(l, b) = 4 and it is easy to verify that the lemma holds. 

• If b = 0, then the induction hypothesis gives 

av'(n, 0) = 1 + av'(n - 1, 1) < 1 + av + (n -1,1) = 

17, . 7 ln-1 17 . 

= H (n - 1) H h 9 = — n + 9 = av + (n, 0). 

12 V J 12 6n- 1 12 v ' ; 

• For b = n we get 

av'(n, n) = 2 + at/(n — 1, n — 1) < 2 + at> + (n — 1, n — 1) = 

17 7 1 17 7 1 

= 2 H (n - 1) H (n - 1) h 9 = — n H n h 9 = av + (n, n). 

12 v ; 12 v ; 6 12 12 6 V ' ; 
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In the case 1 < b < n 

n(n — l)(av + (n, b) — av'(n, b)) 

>n(n — l)av + (n, b) — (n — b)(n — b — 1)(1 + av + (n — 1,6 + 1)) 

- 2(n - b)b ^~ + av + {n - 1, ft) J - b(b - 1) (2 + av + {n - 1, b - 1)) 

6 (\ 1 \ 

-n 6 > 0. 



n - 1 V3 3 , 

□ 

Therefore at> + (n, b) > av'{n,b) and thus 

17 

av'(n, 0) < ai7 + (ra, 0) = — n + 9. 

□ 



6 Computational results 

Computer search found the following sequence of 30 flips that sorts the stack —I\g: (19, 
14, 7, 4, 10, 18, 6, 4, 10, 19, 14, 4, 9, 11, 8, 18, 8, 11, 9, 4, 14, 19, 10, 4, 6, 18, 10, 4, 7, 
14). Thus, using Theorem [TJ #(-/i 9 ) = 30. 

We also computed 120) = 32: From [3j Theorem 7]: g{— 120) < Jig) + 2 = 32. 
From Theorem [TJ g(—I 2 o) > 31 and from Lemma [2] follows that if g(— J20) = 31, then 
each flip of the optimal sorting sequence increases the value of the function v by 4/3. But 
computer search revealed that starting at —I20 we can make a sequence of only at most 
29 such flips. 

The values /(18) = 20 and /(19) = 22 were computed by the method of Kounoike 
et al. [UJ and Asai et al. [TJ. It is an improvement of the method of Heydari and Sud- 
borough [Sj. Let be the set of stacks of n unburnt pancakes requiring m flips to 
sort. For every stack U G U™, 2 flips always suffice to move the largest pancake to the 
bottom of the stack, obtaining stack TJ' . Since then, it never helps to move the largest 
pancake. Therefore TJ' requires exactly the same number of flips as TJ" obtained from TJ' 
by removing the largest pancake and thus TJ" requires at least m — 2 flips. 

To determine for all i G {m,m + 1, . . . , f(n)}, it is thus enough to consider the 
set ^m^=m-2^^n-v ^ n eac ^ stack from this set, we try adding the pancake number n 
to the bottom, flipping the whole stack and trying every possible flip. The candidate 
set composed of the resulting and the intermediate stacks contains all the stacks from 
ufl^UJj. Now it remains to determine the value of f(U) for each stack TJ in the candidate 
set. As in [TTJ and [TJ, this is done using the A* search. 

During the A* search, we need to compute a lower bound on the number of flips needed 
to sort a stack. It is counted differently then in [TTJ and [lj: We try all possible sequences 
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of flips that create an adjacency in every flip. If some such sequence sorts the stack, it 
is optimal and we are done. Otherwise, we obtain a lower bound equal to the number of 
adjacencies that are needed to be made plus 1 (here we count pancake n at the bottom 
of the stack as an adjacency). 

In addition, we also use a heuristic to compute an upper bound. If the upper bound 
is equal to the lower bound they give the exact number of flips. 



n 


m 


m 


n 


m 


m 


n 


m 




14 


13 


30,330,792,508 


15 


15 


310,592,646,490 


16 


17 


756,129,138,051 


14 


14 


20,584,311,501 


15 


16 


45,016,055,055 


16 


18 


4,646,117 


14 


15 


2,824,234,896 


15 


17 


339,220 


17 


19 


65,758,725 


14 


16 


24,974 















Table 2: numbers of stacks of n unburnt pancakes requiring m flips to sort 



Sizes of the computed sets U™ can be found in Table El It was previously known [9], 
that /(18) > 20 and /(19) > 22. No candidate stack of 18 pancakes needed 21 flips thus 
/(18) = 20. Then /(19) = 22 because /(19) < /(18) + 2 = 22. 

The following modification of this method was also used to compute the values of g(n) 
up to n = 17. Again, C™, the set of stacks of n burnt pancakes requiring m flips, is 
determined from the set uffiZm-2^n-v * n a slightly different way. In every stack of 
n burnt pancakes other than —I n (which must be treated separately), some two pancakes 
can be joined in two flips [3, Theorem 1]. We will now show that the two adjacent 
pancakes can be contracted to a single pancake, which decreases the size of the stack. 
The reverse process is again used to determine the stacks of the candidate set, which are 
then processed by the A* search. 

Lemma 10. Let C be a stack of burnt pancakes with a pair (pi,p2) of adjacent pancakes 
and let C be obtained from C by contracting the two adjacent pancakes to a single pancake 
p. Then C can be sorted in exactly the same number of flips as C . 

Proof. If we can sort C in m steps, we can sort C in m steps as well — we do the flips 
below the same pancakes as in an optimal sorting sequence for C. Flips in C below p 
are performed below the lower of pi,p% in C . 

The stack C can be also obtained from C by removing one of the two adjacent 
pancakes. Then we can sort C by doing the flips below the same pancakes as in a 
sorting sequence for C. Flips in C below the removed pancake are performed in C below 
the pancake above it. □ 

During the A* search, we compute two lower bounds and take the larger one. One 
lower bound is computed from the formula in Lemma [3 To compute the other lower 
bound, we try all possible sequences of flips that create an adjacency in all but at most 
two flips. If no such sequence sorts the stack, we obtain a lower bound equal to the 
number of adjacencies that are needed to be made plus 3. 
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In the stacks visited during the A* search, we can contract a block to a single burnt 
pancake thanks to Lemma [T01 If, after the contraction of blocks, the stack has at most 
nine pancakes, we look up the exact number of flips in a table previously computed by a 
breadth- first search starting at Ig. 



n 


m 




n 


m 


1 (P m 1 

l^n 1 


n 


m 


1 (P m 1 

1*^71 1 


n 


m 


|C™| 


10 


15 


22,703,532 


11 


17 


5,928,175 


12 


19 


344,884 


13 


21 


15,675 


10 


16 


179,828 


11 


18 


10,480 


12 


20 


265 


13 


22 


4 


10 


17 


523 


11 


19 


36 


12 


21 


1 


14 


23 


122 


10 


18 


1 














15 


25 


2 



Table 3: numbers of stacks of n burnt pancakes requiring m flips to sort 



Sizes of the computed sets C™ can be found in Table [HJ No stack of 16 pancakes needs 
27 flips thus g(16) = 26 because g(— J 16 ) = 26. Then g(17) = 28 because g(—I 17 ) = 28 
and 0(17) < #(16) + 2 = 28 Theorem 8]. 

The stack obtained from — I n by flipping the topmost pancake is known as J n [3j. Let 
Y n be the stack obtained from — I n by changing the orientation of the second pancake 
from the bottom. The two found stacks of 15 pancakes requiring 25 flips are J15 and 
Y15 and they are the first known counterexamples to the Cohen-Blum conjecture which 
claimed that for every n, —I n requires the largest number of flips among all stacks of n 
pancakes. However, no other J n or Y n with n < 20 is a counterexample to the conjecture. 

Majority of the computations were done on computers of the CESNET METACentrum 
grid. Some of the computations also took place on computers at the Department of 
Applied Mathematics of Charles University in Prague. 

Data and source codes of programs mentioned above can be downloaded from the 
following webpage: http://kam.mff.cuni.cz/~cibulka/pancakes. 

7 Conclusions 

Although the two algorithms presented in Sections H] and have a good guaranteed average 
number of flips, experimental results show that both of them are often outperformed by 
the corresponding algorithms of Gates and Papadimitriou. The average numbers of flips of 
the two new algorithms are very near to their upper bounds calculated in Theorems [6] and M 
and the averages for the algorithms of Gates and Papadimitriou are in Table HI 

We will now design one more polynomial-time algorithm for the burnt version, for 
which no guarantee of the average number of flips will be given, but its experimental 
results are close to the lower bound from Theorem [5j 

Call a sequence of flips, each of which creates an adjacency, a greedy sequence. Note 
that since we are in the burnt version, there is always at most one possible flip that creates 
a new adjacency. In a random stack the probability that we can join the pancake on top 
in a single flip is 50%, therefore starting from a random stack, we can perform a greedy 
sequence of length log 2 n with probability roughly 1/n. The idea of the algorithm is, that 
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whenever we cannot create an adjacency in a single flip, we try all n possible flips and do 
the one that can be followed by the longest greedy sequence. 

As in the previous algorithms, two adjacent pancakes are contracted to a single pan- 
cake. Pancakes 1 and n can create an adjacency (1 is viewed as (n + 1) mod n). Therefore 
when the algorithm obtains the stack (1) we need at most four more flips. 

In Table HI n is the size of a stack, sqp is the average number of flips used by the 
algorithm of Gates and Papadimitriou to sort a randomly generated stack of n unburnt 
pancakes, sgpb is the average number of flips used by the algorithm of Gates and Pa- 
padimitriou for the burnt version and is the average number of flips of the algorithm 
described in this section. 



n 


sgp 


sgpb 


SN 


n + n/ log 2 n 


stacks generated 


10 


11.129 


15.383 


14.935 


13.010 


1000000 


100 


122.925 


150.887 


123.463 


115.051 


100000 


1000 


1240.949 


1502.926 


1127.901 


1100.343 


10000 


10000 


12408.686 


15002.212 


10863.502 


10752.570 


1000 


100000 


124115.000 


150063.000 


106608.900 


106220.600 


10 


1000000 


1241263.600 


1499875.600 


1053866.000 


1050171.666 


5 



Table 4: experimental results of algorithms 



The experimental results together with Theorem support the following conjecture. 

Conjecture 1. The average number of flips of the optimal algorithm for sorting burnt 
pancakes satisfies 

av opt {n) = n + Q 



n 



\ogn 
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